function y = initial_function(x)
if (x >= 0 && x <= pi)
    y = x;
elseif (x >= pi && x <= 2*pi)
    y = 2*pi - x;
else
    y = 0;
end
 

    